Topological defects in two-dimensional crystals 
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By using topological current theory, we study the inner topological structure of the topologi- 
cal defects in two-dimensional (2D) crystal. We find that there are two elementary point defects 
topological current in two-dimensional crystal, one for dislocations and the other for disclinations. 
The topological quantization and evolution of topological defects in two-dimensional crystals are 
discussed. Finally, We compare our theory with Brownian-dynamics simulations in 2D Yukawa 
systems. 
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In 1970's, Kosterlitz and Thouless construct a detailed 
and complete theory of superfluidity on two- dimensions 
(2D) [1]. They indicate vortices pair unbinding will 
lead to a second-order transition in superfluid films. 
Later, a microscopic scenario of 2D melting has been 
posited in the form of the Kosterlitz-Thouless-Halperin- 
Nelson- Young (KTHNY) theory [2]. The KTHNY the- 
ory predicts the unbinding of topological defects to break 
the symmetry in the two-stage transitions. Topologi- 
cal defects, which are a necessary consequence of bro- 
ken continuous symmetry, play a important role in two- 
dimensional phase transition [3, 4]. 

In 2D crystal, the evolution of topological defects have 
been studied in experiments and computer simulations. 
A serial experiments were performed to calculate dis- 
locations and disclinations dynamic of two-dimensional 
colloidal systems, and dissociation of dislocations and 
disclinations were observed [5]. During the years, a large 
number of computer simulations indicated that exist a 
two-stage melting as prescribed by KTHNY theory, how- 
ever, results are still controversial [6]. Our previous work 
found that exist a hexatic-isotropic liquid phase coexis- 
tence during the melting of soft Yukawa systems [7]. By 
Voronoi polygons analysis, the behavior of piont defects 
in the coexistence is very complicated. The evolution of 
topological defects during the melting of two-dimensional 
system still a open question. 

It is interesting to consider the appropriate form for 
the point defect densities when expressed in terms of 
the vector order parameter field <f)(r,i). This has been 
carried out by Halperin [8], and exploited by Liu and 
Mazenko [9]. However, their analysis is incomplete [10]. 
In two-dimensional system, a gauge field-theoretic for- 
malism developed by Kleinert [11]. The gauge theory 
of topological quantum melting in 2+1 dimensions Bose 
system was developed by Nussinov et al, and the Super- 
fluidity and superconductivity can arise in a strict quan- 
tum field-theoretic setting [12]. 

Recently, a topological field theory for topological de- 
fect developed by Duan et al [13]. By using ^-mapping 
method and topological current theory, the evolution of 
topological defect which relate to sigularities of the order- 




FIG. 1: Voronoi polygons analysis in 2D Yukawa systems. 



paramter field, such as vortex in BEC [14] and super- 
conductivity [15], was studied. In this Letter, we will 
develop a topological current theory of discloations and 
disclination in 2D crystals. By using the topological field 
approach, the topological quantization and evolution of 
topological defects in 2D crystals was discussed. We com- 
pared our theory with Brownian-dynamics simulations in 
Yukawa systems. 

In continuum elasticity theory, the homogeneous equa- 
tion in 2D triangular solid is given by [16, 17] 

V 4 x/^ = tkidkdiO + eikd k (ejididjUi) (1) 

where Y = 4/i(/i + A)/(2/i + A) is the two-dimensional 
Young's modulus and <jo = A/(2/i + A) the two- 
dimensional Poisson ratio. The defects associated with 
the continuum elastic theory of a solid are dislocations 
and disclinations. In the following we consider only the 
triangular lattice since it is the most densely packed one 
in two-dimensional and favored by Nature. 

Topological defects associated with the continuum 
elastic theory of a two-dimensional solid are dislocations 
and disclinations (see Fig. 1). Disclinations can be intro- 
duced into the theory in a way similar to the discussion 
of vortices in two-dimensional XY model [18]. Since the 
bond angle is defined only up to 27r/6, it implies 

jd0 = jv0ds = ^-W Q , W 6 = 0,±l,...,n. (2) 
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Where W 6 is the winding number of disclinations. Discli- 
nations, which can influence on the decay of orienta- 
tional order, have a much higher energy than disloca- 
tions. They are defined in terms of the bond angle field 
0(r) = l/2(du y /dx — du x jdy\ which measures the bonds 
orientational order. It is convenient to define an order pa- 
rameter for bond orientations, which for the triangular 
lattices is ip(r) — ipoe 169 ^. However, the bond angle field 
is undefined at the disclination cores, i.e., the zero points 
of the order parameter. We rewrite the orientation order 
parameter ip(r) = + i(f)Q instead of tp(r) = tp^e 169 ^. 
We can define the unit vector field uq (x) as ng = </>g / 1 1 06 1 1 
and 1 1 (f)Q 1 1 = ^^6^6' a = 1? 2. We can construct a topo- 
logical current of the orientation order parameter field 
in two-dimensional crystal, which carries the topological 
information of 4>(x), 
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It is the ^-mapping topological current for disclinations. 
We will see that this current associate with the singulari- 
ties of the unit field n§(x). It is easy to see that this topo- 
logical current is identically conserved, i.e., dkj^ isc = 0. 

We consider only the triangular lattice since it is the 
most densely packed one in two-dimensional and favored 
by Nature. The dislocation is characterized by a Burgers 
vector 6, which can be determined by considering a loop 
enclosing the dislocation. In mathematical language, 
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diUidxi = bi. 
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By introduce two-dimensional local transition order pa- 
rameters as p x = pQx e%GxUx — 4>x + ifix an d Py — 



p<dy e%GyUy — + i<t>yi where G is any reciprocal lat- 
tice vector and po = e lG ° r ° . We can define a unit vector 
field n/(x), where f^(x) = e tGlUl = n\ +inf (I = x,y). It 
can be further expressed as 
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||^(x)|| = ^f(x)^(x), 



obviously, nfnf = 1. We can construct topological cur- 
rents of the order parameter field in two-dimensional 
crystal, which carries the topological information of 
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where the Burgers vectors satisfy b\G\ = 2tt. It is the 
^-mapping current for dislocations. Dislocations can be 
quite effective at breaking up translational order. How- 
ever, they are less disruptive of orientational correlations. 
The disclination can be seen as a bound of dislocations 
with parallel Burger vectors. On the hand, disclinations, 
which play a important role in the disruptive of the ori- 
entational order in two-dimensional solid, can be seen as 
bound disclination and anti-disclination pairs. 



From the relationship between the bond angle field and 
displacement vector field that diO = —e a bd a di(u b ). Using 
the topological current theory, we can obtain 



dkjdlso 



~ € HJdi. 



(6) 



it indicate that disclinations can be seen as sources for 
dislocations. Since there aren't free disclinations in two- 
dimensional crystal at the low temperature phase, the 
dislocation currents are conserved as dk J\. = 0. 

By using Green function relation and substitute 0g to 
the homogeneous equation [13], we can rewrite the ho- 
mogeneous equation into a (5-function form as 

>x = ^ 2 WJ(^) + e «4^ 2 («J(^). (7) 

Y O X Lrl X 

With the ^-function theory, the homogeneous equation 
can deduce into 
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(8) 

where r\i = ±1 is Brouwer degree, the positive integer $ 
is called the Hopf index of map x — > (j). The topologi- 
cal charge of disclination and dislication is fyrji = 2tt/6 
and /3jT]j = b l y The first term is related to dislocation, 
and the second term is related to disclination. Also, it 
can simply seen as disclination current, if dislocations are 
thought of as disclination dipole pairs. 

We explore what will happen to the disclination at the 
zero point According to the values of the vector 

Jacobian at zero points of the order parameter, there 
are limit points (J°((/)q/x) ^ 0) and bifurcation points 
(J°((/)e/x) = 0). Each kind corresponds to different cases 
of branch processes. The limit points are determined by 
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Considering the condition (9) and making use of the im- 
plicit function theorem, the solution in the neighborhood 
of the point is t = t(x), y = y(x), where t* = 

t(zj). In this case, one can see that dx/dt\^ ^ = oo, 
the Taylor expansion of t = t(x) at the limit points (£*, ij) 
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(10) 



which is a parabola in the x-t plane. From this equa- 
tion, we can obtain two solutions x\(t) and ^(t), which 
give two branch solutions (World lines of disclinations). 
If d 2 t/ dx 2 \( t * ^ > 0, we have the branch solutions for 
t > t* . It is related to the origin of a disclination 
pair. Otherwise, we have the branch solutions for t < £*, 
which related to the annihilation of a disclination pair. 
Since the topological current is identically conserved, the 
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topological charges of these two generated or annihilated 
disclinations must be opposite at the limit points, say 
PiVi + P2V2 = 0. For a limit point it is required that 
J 1 (^)lf,ii^0. 

At the neighborhood of the limit point, one can obtain 
the scaling law as [19] v oc (t — t*) -1 / 2 . The growth 
or annihilation is parameterized in terms of a relevant 
characteristic length £(£), which is also characteristic the 
mean separation distance of disclinations pair. The rel- 
evant length £(t) obeys that ~ (t - T) 1 / 2 . It is 
the dynamic scaling law of the topological defects pairs. 
In the low temperature equilibrium phase has essentially 
no disclination pairs, and £(£) is infinite. The relation- 
ship between the relevant length and the metastable 
disclination density which below the critical temperature 
is£(t) = y/l/py. then the number of disclinations satisfy 
the power law N oc t -1 . 

Consider in which the restrictions on zero point (£*, z\) 
are J k ((j)§/x)\^ ^ = 0, {k = 0, 1, 2), which imply an im- 
portant fact that the function relationship between t and 
x or y is not unique in the neighborhood of the bifurca- 
tion point (t*,zi). With the aim of finding the different 
directions of all branch curves at the bifurcation point, we 
suppose d(j)\/dy\ t ^ ^ 7^ 0. According to the ^-mapping 
theory, the Taylor expansion of the solution of the ze- 
ros of the order parameter field in the neighborhood of 
(£*, z{) can be expressed as 
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A(x-z}) 2 + 2B(x-zl)(t-t*)+C{t-t* 



0, (11) 



where A, B, and C are constants determined by the order 
parameter. The solutions of Eq. (11) give different direc- 
tions of the branch curves (world line of disclinations) 
at the bifurcation point. There are four possible cases, 
which will show the physical meanings of the bifurcation 
points. 

First we consider the case that A ^ 0. If A = 
4(B 2 —AC) > 0, from Eq.(ll) we get two different motion 
directions of the core of disclination dx/dt\ 12 = (— B ± 

\J B 2 — AC) /A, where two topological defects encounter 
and the depart at the bifurcation point. However, when 
A = 4(B 2 - AC) = 0, form Eq. (11), we obtain only one 
motion direction of the core of disclinationcb /d£ | 2 = 
—B/A, which includes three important cases, (i) Two 
disclinations tangentially encounter at the bifurcation 
point, (ii) Two disclinations merge into one disclination 
at the bifurcation point, (iii) One disclinations splits into 
two disclinations at the bifurcation point. 

Then we consider the case that A = 0, If C 7^ 
and A = 4(£ 2 - AC) = , from Eq. (11) we have 
dt/dx\ 12 = (-B ± VB 2 - AC)/C = 0, -(2B)/C. There 
are two important cases, (i) One disclination split into 
three disclinations at the bifurcation point, (ii) Three 
disclinations merge into one disclination at the bifurca- 
tion point. From the above discuss, we can then ob- 
tain the growth or annihilation velocity of the vortices as 
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FIG. 2: (a)The fractions of 6-coordinated, 5-coordinated, and 
7-coordinated particles. At low temperature the number of 
defects is very small. When the temperature above 0.500, 
N&/N decreases rapidly. 



v oc const. The approximation asymptotic relation of is 
£(t) oc (t - T). If C = 0, Equations (11) give dx/dt = 
or dt/dx = 0. This case shows that two worldlines inter- 
sect normally at the bifurcation point, which is similar 
to the case that A = and C / 0. We can obtain 

= const, v = 0. 
It is obvious that vortices are relatively at rest when 

= const. Through our topological current theory, we 
can get the dynamical scaling law. Moreover, it worth to 
emphasizing that the dynamical scaling law of vortices, 
only depend on topological properties of the order pa- 
rameter field. 

We performed Brownian dynamics simulations on 
melting of two-dimensional colloidal crystal in which par- 
ticles interact with Yukawa potential [7]. We have ob- 
served the mechanisms of defects in the two-dimensional 
Yukawa system [20]. Fig. 2 plots the fractions of 
6-coordinated, 5-coordinated, and 7-coordinated parti- 
cles. At low temperature, all particles are nearly six- 
coordinated, and the number of defects is very small. 
When the temperature reaches 0.500, Nq/N behaves in 
a rapidly decreasing. At T* = 6.05, almost 20% of parti- 
cles are attached to the defects, which consisted with the 
results in Ref. [? ]. As the defects fraction rises above 
30%, the system melt into a liquid phase. At the low 
temperature, the paired dislocation is formed or annihi- 
lated. With the temperature rises to the hexatic phase, 
the dislocation pair comes into dissociated. In the solid 
phase, the unbinding of dislocation is unstable, and these 
unbind dislocation will bind soon, whereas the stable free 
dislocation is found in the hexatic phase. The above solu- 
tions form topological current theory reveal the evolution 
of disclinations. 



4 



Besides the encountering of disclinations, splitting and 
merging of disclinations are also included. When a mul- 
ticharged disclinations, such as 8-coordinated particles, 
pass the bifurcation point, it may split into several discli- 
nations with lower value of Burgers vector along different 
branch curves. On the other hand, several disclinations 
can merge into one disclination at the bifurcation point. 
As before, since the topological current of disclinations is 
identically conserved, the sum of the topological charge 
of final disclinations must be equal to that the initial vor- 
tices at the bifurcation point, i.e. ^fPifrjif = ^iPurju 
for fixed 1. It indicates that, vortices with a higher value 
of Burgers vector can evolve to the lower value of Burgers 
vector, or vortices with a lower value of Burgers vector 
can evolve to the higher value of Burgers vector through 
the bifurcation process. But due to Such case was ob- 
served in our BD simulation in Yukawa systems. That 
is why at the liquid phase the number of 5-coordinated 
particles is much more than 7-coordinated particles due 
to the emergence of the 8-coordinated particles. The sim- 
ilar experimental result was observed in melting of two- 
dimensional tunable-diameter colloidal crystals [21]. 

We would like to thank Y. X. Liu for his helpful discus- 
sions. This work was supported by the National Natural 
Science Foundation of China and by the SRF for ROCS, 
SEM. 



* Email: ychen@gmail.com 

* Email: weikaiqi@gmail.com 

[1] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973). 

[2] D. R. Nelson and B. I. Halperin, Phys. Rev. B. 19, 2457 
(1979); A. P. Young, Phys. Rev. B 19, 1855 (1979); K. 
J. Strandburg, Rev. Mod. Phys. 60, 161 (1988). 

[3] D. R. Nelson, Defects and Geometry in Condensed Mat- 
ter Physics, (Cambridge University Press, Cambridge, 
2002). 

[4] P. M. Chaikin and T. C. Lubensky, Principles of con- 
densed matter physics, (Cambridge University Press, 



Cambridge, 1995). 
[5] C. A. Murry and D. H. Van Winkle, Phys. Rev. Lett. 58, 
1200 (1987); Y. Tang, A. J. Armstrong, R. C. Mockler, 
and W. J. O'Sullivan, Phys. Rev. Lett. 62, 2401 (1989); 
A. H. Marcus and S. A. Rice, Phys. Rev. E 55, 637 

(1997) ; H. H. von Griinberg, P. Keim, and G. Maret, 
Soft Matter ( Vol.3): Colloidal Order from Entropic and 
Surface Forces, Edited by G. Gompper and M. Schick, 
Wiley- VCH (2007). 

[6] K. Chen, T. Kaplan, and N. A. Clark, Phys. Rev. Lett. 

74, 4019 (1995); K. J. Naidoo and J. Schnitker, J. Chem. 

Phys. 100, 3114 (1994). 
[7] W. K. Qi, S. M. Qin, X. Y. Zhao, and Y. Chen, J. Phys.: 

condensed Matter 20, 245102 (2008). 
[8] B. I. Halperin, Physics of Defects, edited by R. Balian et 

al. (North-Holland, Amsterdam, 1981). 
[9] F. Liu and G. F. Mazenko, Phys. Rev. B 46, 5963 (1992). 
[10] Y. S. Duan, H. Zhang, and L. B. Fu, Phys. Rev. E 59, 

528 (1999). 

[11] H. Kleinert, Gauge Theory in Condensed Matter, (World 

Scientific, Singapore, 1989). 
[12] J. Zaanen, Z. Nussinov and S. I. Mukhin, Ann. Phys 

(New York) 310, 181 (2004); V. Cvetkovic, Z. Nussinov, 

S. Mukhin, and J. Zaanen, Europhys. Lett. 81, 27001 

(2007). 

[13] Y. S. Duan and S. L. Zhang, Int. J. Eng. Sci. 29, 689 
(1991); G. H. Yang and Y. S. Duan, Int. J. Theor. Phys. 
37, 2371 (1998); Y. S. Duan and H. Zhang, Phys. Rev. 
E 60, 2568 (1999). 

[14] Y. S. Duan and H. Zhang, Eur. Phys. J. D 5, 47 (1999). 

[15] Y. S. Duan, H. Zhang and S. Li, Phys. Rev. B. 58, 125 

(1998) . 

[16] L. D. Landau and E. M. Lifshitz, Theory of Elasticity, 

(Pergamon, New York, 1970) 
[17] H. S. Seung and D. R. Nelson, Phys. Rev. A 38 1005 

(1984). 

[18] W. K. Qi and Y. Chen, e-prints arXiv: 0809.0348. 

[19] T. Xu, Phys. Rev. E. 72, 036303 (2005). 

[20] A movie is archived at the website: http : 

// '210. 26. 51. 253 /tiki /tiki — download file, phpl f Held — 

17. 

[21] Y. Han, N. Y. Ha, A. M. Alsayed, and A. G. Yodh, Phys. 
Rev. E 77, 041406 (2008). 



